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Abstract 

Density-dependent relations among saturation properties of symmetric nuclear matter and 
hyperonic matter, properties of hadron- (strange) quark hybrid stars are discussed by applying 
the conserving nonlinear a-u-p hadronic mean-field theory. Nonlinear interactions that will be 
renormalized as effective coupling constants, effective masses and sources of equations of motion 
are constructed self-consistently by maintaining thermodynamic consistency to the mean-field 
approximation. The coupling constants expected from the hadronic mean-field model and SU(6) 
quark model for the vector coupling constants are compared; the coupling constants exhibit dif- 
ferent density-dependent results for effective masses and binding energies of hyperons, properties 
of hadron and hadron-quark stars. The nonlinear a-uj-p hadronic mean-field approximation with 
or without vacuum fluctuation corrections and strange quark matter defined by MIT-bag model 
are employed to examine properties of hadron- (strange) quark hybrid stars. We have found 
that hadron- (strange) quark hybrid stars become more stable in high density compared to pure 
hadronic and strange quark stars. 

PACS numbers: 21.65. +f, 24.10. Cn, 24.10. Jv, 26.60.+C 

1 Introduction 

The symmetric nuclear matter is a self-bound matter whose binding energy exhibits a charac- 
teristic concave curve at saturation density, where pressure of nuclear matter vanishes (p = 0). 
It has been known as a constraint to examine self-consistency to nuclear many-body approxi- 
mations pQ — [3]. The self-consistency has been fundamental for many-body approximations in 
terms of Landau's Fermi-liquid theory [3] — [B], Kadanoff-Baym's theory of conserving approxima- 
tions [7J — [9] and the density functional theory |10j — [12]. The relativistic, field-theoretical, linear 
a-oj mean-field approximations that maintain conditions of conserving approximations have been 
applied to finite nuclei, nuclear matter and neutron stars \Y'i\ I14| . These nuclear mean-field ap- 
proximations reveal an important relation among self-consistent single particle energy, energy 
density and particle density, which is denoted as thermodynamic consistency |15t 116] . The con- 
serving nonlinear a-uo-p mean-field approximation which maintains thermodynamic consistency 
has been applied to investigate density-dependent correlations among properties of nuclear, 
neutron and hyperonic matter [T7] — [TO]. Thermodynamic consistency is important to derive 
consistent results for high energy and high density phenomena, such as heavy-ion collision, 
hypernuclei- formation experiments [20] and astrophysical phenomena [21] — [23]. 

The nonlinear a-uo-p mean-field approximation has several nonlinear coefficients whose values 
are determined by reproducing binding energy, symmetry energy and simultaneously searching 
the minimum value of incompressibility at saturation density, as well as reproducing a maximum 
mass of neutron stars at high density. The binding energy at saturation is taken as —15.75 MeV 
at k F = 1.30 fm _1 (po = 0.148 fm~ 3 ) with the symmetry energy, 04 = 30.0 MeV in the current 
calculation |13] , The values of nonlinear coefficients which will produce effective masses of 
hadrons M* N /M N ~ 0.70, m* a /m a ~ 1.02, m* /m, ~ 1.01 and incompressibility, K ~ 320 MeV, 
with symmetry energy 04 ~ 30 MeV, the maximum mass of neutron stars in /3-equilibrium 
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matter [23], M max (n,p,e) ~ 2.50 M (the solar mass: M Q ~ 1.989 x 10 30 kg), are found 
reasonable and admissible in the nonlinear a-uo-p Hartree approximation |181 I19j . 

The maximum masses of neutron stars are expected to be below 2.5 M [24] . or 2.1 ± 
0.2Mq |25j . which depends on the hadronic equation of state (EOS) for isospin asymmetric, 
hyperonic matter in the density range, 2pQ < p B < 5po. The hadronic liquid-gas phase transition 
at the surface of neutron stars and determination of the radius of smeared, gas-phase surface are 
not directly relevant to properties of hadronic matter and maximum masses of neutron stars. 
The nonlinear interactions exhibit significant density-dependent effects on incompressibility, K, 
and symmetry energy, 04, in high densities; these Fermi liquid properties monotonically increase 
about saturation density, but they are piecewise continuously softened at an hyperon onset 
density. These phenomena of piecewise continuous change of Fermi-liquid properties will be 
important for the analysis of Landau parameters, heavy-ion collision, high energy and high 
density experiments. At a hyperon onset density from (n,p, e) to (n,p, H, e), the EOS suddenly 
becomes softer. This is because the nucleon Fermi energy, E^(k F ) in the phase (n,p, e), will 
be redistributed to the hyperon Fermi energy, En{k F ) in the phase (n,p,H,e); consequently, 
the Fermi energies become relatively small in the phase (n,p, H,e) compared with those of 
(n,p,e). The redistribution and slow increase of Fermi energies appear whenever hyperons are 
generated, resulting in a softer equation of state and discontinuous changes of K and 04. This 
is numerically checked as discrete changes of physical quantities, such as effective masses of 
hadrons, incompressibility, symmetry energy and energy density |X8|, I19j . 

Since the hyperon-onset will confine Fermi-energies of baryons as explained above, single 
and double hyperon generations exhibit different density-dependent phenomena. For example, 
the A-hyperon onset density such as in (n,p, e)-(n,p, A, e) is p\/po ~ 2.2; however, it gives 
Pa/ Po ~ 4.2 when A is generated in (n,p, e)-(n,p, E~, e)-(n,p, X - , A). The same phenomena 
are observed with other hyperons, and generally the onset-density of a hyperon is pushed up 
to a higher density, which is denoted as the push-up phenomena of hyperon onset-densities in 
many- fold hyperon generations [18]. Because of the push-up phenomena of the hyperon onset 
density, we have found that hyperon generations are suppressed in high densities and hyperons 
relevant to determine the maximum mass of neutron stars are E~ and A. The similar results are 
discussed in nonrelativistic Brueckner-Hartree-Fock [261 [27] and quark matter calculations [28] . 
The hyperon onset densities are related to hyperon single particle energies by way of phase equi- 
librium conditions. Hence, the Hugenholtz-Van Hove theorem and thermodynamic consistency 
at saturation of hyperon binding energy are essential to define self-consistent approximation of 
hyperonic matter. The single particle en erg ies of h yperons are also important to study K and 
04 for magic nuclei [29] and hypernuclei [30] — [32] . 

Density-dependent effective masses and effective coupling constants, saturation properties for 
nucleons and hyperons are discussed in the nonlinear a-uo-p mean-field approximation |18[ 119] . 
The results show; (1) coupling constants of hyperons are related to those of nucleons by effective 
masses, effective coupling constants, binding energy at onset-density of respective hyperons. 
Hence, it shows that binding energies of symmetric nuclear matter and hyperon matter are 
self-consistently related to each other. (2) Self-consistency suppresses hyperon generations in 
high densities, denoted as the push-up phenomena of hyperon onset densities. The suppression 
of hyperon generations is also discussed in different calculations [26] — [28]. (3) Coupling ratios 
of hyperons are expected to be g aH j 'g a N,9uiH 1 1 9ujn ^ 1> in order to be consistent with conditions 
of thermodynamic consistency, empirical values of nuclear matter and neutron stars. In the 
current calculations, we have included vacuum fluctuation corrections (VFC) into the nonlinear 
a-u-p approximation and examined properties of (n,p, e), (n,p, H\,e), (n,p, Hi, H2, e) hyperonic 
matter, constrained by transitions to strange quark matter and properties of hadron-quark 
hybrid stars. The hadron-quark phase transitions are assumed to be a first order and computed 
by Maxwell construction [33] . 

The conserving nonlinear a-uo-p mean- field approximation [T7J — [19] and quark-based effec- 
tive nuclear models [M] — [36J have been applied to finite nuclei, nuclear and isospin asymmetric, 
high-density matter. We have compared hyperon coupling constants required by the nonlinear 
a-uj-p mean-field model with those required by SU(6) quark model for the vector coupling con- 
stants [37] [38]. The hyperon coupling constants required by hadronic and SU(6) quark models 
exhibit quite different results for effective masses, binding energies of hyperons [181 02] an d 
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properties of hadron-quark hybrid stars. The coupling ratios required by SU(6) quark model 
cannot reproduce hyperon saturation properties, which will be discussed in terms of conditions 
of thermodynamic consistency in the sec. 6. The analysis of discrepancies of predictions by 
hadronic and quark-based models may provide us with insight for constructing self-consistent 
nuclear many-body problems [UJ — |19j . The effective masses of hyperons depend on coupling 
ratios (g aH / 1 9oN = r^ N , g u H / 9uN = r%^, g P H/g P N = r p HN ) and decrease analogous to effective 
masses of nucleons in high densities, which shows strong density-dependency of hadronic inter- 
actions. However, the effective masses of hyperons exhibit weak density-dependent results with 
r HN — 2/3 required by SU(6) quark model for the vector coupling constants. In general, hadrons 
exhibit strong density-dependent interactions and correlations among properties of nuclear mat- 
ter, hyperonic matter and neutron stars. In the hadronic mean-field approximation, rfj N > 1.0 
is preferred in order to be consistent with properties of nuclear matter and the maximum mass 
of isospin asymmetric neutron stars, M max (n,p, e) = 2.00 ~ 2.50 M . The value, r L f IN = 2/3, 
results in different effective masses and binding energies of hyperons; the discrepancies originate 
from density-dependent interactions of hadrons [19J . 

The density-dependent many-body effects produced by the conserving nonlinear a-u-p mean- 
field approximation should be compared to chiral hadronic mean-filed approximations [35] — |42j . 
The effective quark-based chiral lagrangian approach suggests that the appropriate in-medium 
scaling law, m*/m a « m* N /miy ~ m^/m^ ss /*/ 'f w m*/m p , be expected [35j 03]. Although 
it is ambiguous whether m* N /mN m m£,/m w means that m*/m u < 1.0 or m*/m w > 1.0 at 
saturation density, they are certainly decreasing above saturation density, since m* N /rriN is 
model-independently expected to be decreasing. This is significantly different from nonlinear 
mean-field approximations, since self-consistency of hadronic mean-field approximation demands 
m^/m a w m^j/rriuj > 1.0 [T71 05] and m*/m p rn^/m^ < 1.0 [36] in high densities. In 
other words, since effective masses and coupling constants are self-consistently related to each 
other, if some values decrease, the others have to counterbalance the variations. Therefore, 
all correlated effective masses and coupling constants must decrease or increase simultaneously 
in order to compensate for variations among others. At saturation density, effective masses of 
mesons are density-dependent and their ratios are approximately equal: m*/rra CT ~ w.*/m w ~ 
m*p/m p ~ m^/m-x ~ 1.0. However, mesons exhibit different behavior above saturation densities 
as explained. The discrepancy has been discussed in terms of thermodynamic consistency in 
the hadronic mean-field approximations [17] — [TU] and should be investigated in other hadronic 
models to extract consistent model-independent results. 

We have applied the nonlinear a-u-p mean-field approximation and MIT-bag model upon 
hadron-quark hybrid stars [231 07] . The numerical analysis exhibits new results that the hadron- 
(strange) quark hybrid stars are more stable in high density than pure hadronic and strange 
quark stars. It suggests a relation between bag constant and QCD strong coupling constant, 
(B, a c ) [481 149] . to the central density and maximum mass of hadron-(strange) quark hybrid 
stars, (£ c , M max ). The results obtained in the current calculations should be examined from 
astronomical data whether or not hadron-quark stars are possible and the values of bag constant 
and strong coupling constant, (B, a c ), could be consistent with astronomical data for neutron 
stars. 

Self-consistent relations among saturation properties of nucleons and hyperons are briefly 
reviewed in sec. 2. Quantitative numerical calculations for effective masses, onset densities and 
conditions of hyperon saturation are discussed in the articles [181 IS]- The MIT-bag quark 
matter and vacuum fluctuation correction to nonlinear mean-field approximation are explained 
in sec. 3. Results of pure hadron and hadron-strange quark hybrid stars are discussed in sec. 4; 
concluding remarks are in sec. 5. 

2 Self-consistent effective masses and coupling constants in the 
nonlinear a-tu-p mean-field approximation 

The hadronic lagrangian with nonlinear a-u-p in tera ctions which yields density-dependent ef- 
fective masses and coupling constants is given by [17J, 
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—L^ ■ IT + l -m 2 p R^ ■ R» + 9 -^-{R„ • i^) 2 + ^ 2 i^ • + • R? 

+ ^i>i{il^-m l )^ l + 5C (2.1) 
i 

where ips (B = n,p, A, S, •••) and t/>; (Z = e~ , fi~) denote the field of baryons and leptons, 
respectively. The meson-fields operators are: <j) for the cr-field, V for the vector-isoscalar oo- 
meson, V^V^ = Vq — V 2 , (/x = 0,1,2,3) and R p for p-meson. The vector field strengths, F pv 
and L pu , are defined as, F pu = d p V u - d u V p and L pu = d^Ry - d v R p + g p Rfj,x R v . 

The coupled nonlinear quantum-field lagrangian (|2.ip is interpreted as baryon quantum-field 
lagrangian in mean fields of mesons [13]. All dynamics generated by baryon fields are mediated 
by mean fields of mesons which will be self-consistently defined in an approximation. The ef- 
fective coupling constants, g* B (i = a,u:,p), denote renormalized, density-dependent coupling 
constants defined by self-consistent mean-field of u-meson. The nonlinear a-uj-p mean-field 
model maintains the structure of Serot and Walecka's linear a-uj mean-field approximation [13] . 
Lorentz-invariance and renormalizability, the Hugenholtz-Van Hove theorem [2j, conditions of 
conserving approximations [7J — [12], the virial theorem |151 [T6] and Landau's hypothesis of 
quasiparticles [1] — [6]. As we proved in the ref. [17], nonlinear mean- field approximations are 
equivalent to Hartree approximation when nonlinear interactions are properly renormalized, 
and consequently, the concepts of effective masses and effective coupling constants are naturally 
generated by nonlinear mean- field interactions. Self-consistent relations among single parti- 
cle energy, effective masses and coupling constants will restrict empirical values of low-density 
nuclear matter and high-density hadronic matter. The admissible values of effective coupling 
constants and masses are confined in certain values due to strong density-dependent correlations 
among physical quantities of nuclear matter and neutron stars [T7J — |19j . 

Meson-fields operators are replaced by mean- fields denoted by 4>q, Vq, and Rq- The equations 
of motion for baryons are given by, 

(i 7 ^ - g^ioVo - ^%T 3 #o) - (M B - g* aB cf> )] ips = , (2.2) 

where g* B , g^ B and g* B are effective coupling constants for a, uj and p mesons. One should 
notice that effective coupling constants cannot be simply introduced as experimental, external 
inputs to an employed approximation, since density-dependent coupling constants will modify 
equations of motion for mesons. We will assume that only nucleon-meson coupling constants are 
density-dependent since we are interested in the density-dependent correlations among properties 
of symmetric nuclear matter. Nonlinear interactions are not assumed in the coupling constants 
for hyperons; effective masses of hyperons are defined by Mjj = Mh — g<jH<fio- The density- 
dependent nucleon-meson coupling constants that maintain thermodynamic consistency are, 

9a N = 9a N + gaaN(f>o/2rn a , 

9wN = 9wN + 9aLjN4>o/m a , (2.3) 
9* p n/ 2 = 9 p n/1 + ga P N<t>o/m a . 
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The equations of motion for mesons are given by, 

2 i , 9o3 ,2 , 9oA ,3 9uui^ r 2i 9a p n 2 i 9aaN , * 9cruN lr 9apN n 

m^Po + -y-(f> + -gj-^o - ^-^o 'Po - -^o^o - 2^ _/9s ^° = 5ctJvPs ~ "T^ - oAj ~ 0p3 ' 

m 2 p Ro + 9 fRl + 9 ~f4Ro + 9 -fViR = 9 -f P s . (2.4) 

where p s = J2b PsB is the total scalar source, p u the isoscalar density, and p 3 = (k Fp — k Fn ) /3tt 2 , 
the isovector density. The density-dependent coupling constants will modify the equation of mo- 
tion for u-meson which acquires a mass term, 9 £££ p s 4>o, and new source terms, —gauNVoPu/m^ — 

g CTp NRoP3/nT'a from density-dependency of effective coupling constants. 

The introduction of nonlinear UCT./V- vertex interaction leads to the effective mass of nucleoli: 

Ml; = M N - g* aN (po = M N - swv^o - (flwjv /2m a )4>l > ( 2 - 5 ) 
and effective masses of nucleons and hyperons are related to each other as, 

M H -M* H = ?^(M N -M* N ) ■ (2.6) 

9oN 

The total scalar source is obtained by the requirement of self-consistency, 

S s = Z% + E^ = -%(p* sN + Psh) , (2.7) 
where the scalar sources are respectively given by 

g* B /g* N [ kF *i2 M* B 



Jo E* B {q) 

and p* N is the modified scalar density defined by g* N p* sN = g* N PsN-gaujNVoPB/m a -g apN RoP3/m a ; 
N is used to denote proton and neutron, N = (p, n); hyperons are denoted as, H = A, E~, E°, E + , • • • 
The w-meson and />meson contributions to self-energies are given by 

^ = -ffM,o and ^pT^, (2-9) 
where the isoscalar density, p w , is given by 

Pu = Pp + Pn + ^YjTHNPH, (2-10) 
H 

and r BN = g^H I 'g^N 1S ^ ne density-dependent ratio of hyperon-nucleon coupling constants. 
The self-energies, Ep„ and Ep n , are briefly denoted as E, p ,; the isovector density is denoted as 

P{n) 

P3 = (kp — kp n )/3ir 2 where the Fermi momentum kp p is for proton and fcp n for neutron |17t I18J . 

The baryon-isovector density, p 3 B, and the ratios of sigma-nucleon coupling constants on p- meson 
are also defined; for example, p 3B = p 3 + r^ N p 3j: , where r^ N = g p z/g* N and /9 3E = p s + - p s -. 
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Thermo dynamically consistent effective masses of mesons compatible with effective coupling 
constants (|2.3p are required to be: 



in 



*2 



2 / r i So-3 
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*2 



m 
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2m 2 
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+ 
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3!m 



2^0 
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m 2 / i , flp4 „2 , ffgp .2 , T/ 2 

mj 1 + tti-^-Ko + TTZ^O + 2^2 
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2m 2 



gap 

2m 2 



9auN 
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■PsN 



(2.11) 



Since effective masses of mesons and coupling constants depend on mean fields of mesons, they 
are density-dependent through meson fields. Note that the effective mass of a-meson depends on 
the (n,p) scalar source of nucleons, p s rv- The modifications to equations of motion, propagators 
and self-energies produced by density-dependent effective coupling constants and masses have 
to be carefully discussed. 

The energy density, pressure of isospin-symmetric, asymmetric and charge- neutral hadronic 
matter are calculated by way of the energy-momentum tensor: 



B 



Tr- 
ill 
~2 



kFB ,i 1 2 T p n\ i m o- j2 , 9cr3 , 3 . 9crA , 4 m t^ r 2 9u4 17 -4 9aw 

dkk E B (k) + —<j) Q + -y-(f> Q + —(p Q - —V - —V - —c 



2 , 9crp 

o + -r ( 



E 



1 



7T^ 



dkk 2 Ei{k) , 



(2.12) 



Pnha 



3tt 2 



dA; 



A' 1 



E* B (k) 



m l A2 S<t3 ^ 3 9*4 ^ ™l 2 



3! 



4! -o^— "o ^Tf^o 



F 4 , 5^12^2 



2 

/ m p 9p4 R 2 , 9ap ,2 9u)p T/ 2\ p2 , 1 f 



^ „ A; 4 
dk 



l=e ,p 



E*{k) ' 



(2.13) 



where kF B is the Fermi-momentum; Esik) and Ei(k) are single particle energies for baryons and 
leptons, respectively. One can check that the thermodynamic relations, such as £mha -\~Pnha — 
p B E n (kF n ) and the chemical potential, p = d£ NHA /dp B = E n (kF n ) = E*(kp n ) — T, (kp n ), are 
exactly satisfied with a given baryon density, p B = 2kp/3ir 2 . 

The functional derivative of energy density, 8 NHA ((f>o, Vq, Rq, th), with respect to the baryon 
number distribution, nj, is given by: 



88, 



5m 



E(h)+J2 



58, 



5<fto 5n 



o 68 



Nil A 

6V 58^ 



5_Rq \ 



8Vo 8rii 5Rq Srii ) 



(2.14) 



Thermodynamic consistency requires: 



S8 A 



0. 



88„ 



and 



88 a 



[H]. Now, 



acpo 8Vq 8R 

one can directly prove that the self-energies calculated by propagators and the conditions of 
conserving approximations become equivalent, only if the effective masses and effective coupling 
constants of mesons are given by (|2.3p and (|2 lip [17] — [19]. 

In order to start self-consistent hadronic matter calculations, nonlinear coupling constants, 
9aNi 9uiNi 9pN and other 9 nonlinear coefficients should be supplied; but admissible values 
of nonlinear coefficients are determined to satisfy properties of symmetric nuclear matter at 
saturation (—15.75 MeV at k F = 1.30 fin -1 , 04 = 30.0 MeV), simultaneously searching for the 
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minimum value of incompressibility and reproducing the maximum mass of isospin-asymmetric 
neutron stars (M max (n,p, e) = 2.50 M & ). In addition, hadronic phase transitions from (n,p, e) to 
(n,p, H\,H2, ■ ■ ■ , e) should be carefully incorporated in EOS. The binding energies of hyperons 
required by hadronic model and SU(6) quark model for vector coupling constants [371 |3B] are 
compared and discussed quantitatively in the articles [181 fl9] . 

3 The nonlinear a-u-p mean-field approximation with vacuum 
fluctuation correction (VFC) 

We have included vacuum fluctuation correction (VFC) into the nonlinear a-oj-p mean- field ap- 
proximation and applied the EOS with or without VFC to hadronic stars and hadron-quark 
hybrid stars. The EOSs for hadrons with VFC and quark matter generated by MIT-bag model 
are briefly discussed; based on the formalism, the pure hadronic and quark stars, stability of 
hadron-quark stars, are discussed in sec. 6. The vacuum fluctuation corrections are explic- 
itly performed with counterterms required by power-counting and the method of dimensional 
regularization |13j . 

The self-energies with VFC in the conserving nonlinear a-u-p approximation are given with 

by, 

Efj =M* B {k F ) - M B 

1 r / M* \ 5 

=X% + Z S H + [Mf In (jgj ~ M 2 B {M* B - M B ) - -M B (M% - M B f 



^(M* B -M B f] 



(3.1) 



where B = n,p, A, S and g* B = g^ B for B = A, E , since we are investigating 

effects of density-dependent interactions of nucleons; the density-dependent nonlinear self and 
mixing interactions among hyperons will be studied in the future. The self-energies E 1 " and E° 
obtain no VFC in the mean-field approximation, which will be proven directly by dimensional 
regularization method. Note that it is also essential for an a pprox imation with VFC to maintain 
conditions of conserving approximations in order to obtain (|3.ip . 
The energy density with VFC is derived as, 



1 r / M* \ 7 

A£ vfc = -^E [ M b ln + m b(M b - M B ) - —M B (M B - M5) 2 



1 3 95 
+ -M B (M B - M* B f - -(M B - M%f 



(3.2) 



Then, the total energy density is given with eq. (|2,12p as, 

£vfc = £nha + A£ VFC ■ (3.3) 

The self-energies and pressure are also evaluated by dimensional regularization and thermody- 
namic consistency can be proved including VFC. As discussed in sec. 2, all coupling constants 
have to be evaluated by maintaining properties at nuclear matter saturation, searching the min- 
imum value of incompressibility. Incomressibility is K ~ 350 MeV and M max = 2.33 Mq for 
(n,p, e)+VFC matter. The vacuum fluctuation corrections to nonlinear a-uj-p mean-field ap- 
proximation are not significant in low and high densities compared with those of nonlinear a-u-p 
self and mixing interactions. 
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The energy density and pressure of MIT-bag model are derived [17] — |49| : 



- m}k f E f (k f ) - mjlog \ kf + 



+ B, 



(3.4) 



JL £ [ (2 *J _ 3k f m})E f (k f ) + 3m) log 



kf + E f (k f ) 
m f 



B 



(3.5) 



where Ef(kf) = (k'j + m 2 ) 1 / 2 , (f = u,d,s), and I? is the bag-constant (MeV/fm 3 ); fc/ is the 
Fermi-momentum of flavor /. The energy density and pressure with the following given baryon 
density, charge neutrality and phase equilibrium conditions: 



Pd 



214 

3vr 2 



1 



(flu + Pd + Ps 



= 2p u - p d - 
E d (k d ) = E s {h 

determine (u, d, s)-quark matter uniquely (note: pf — n,^ 
for / to first order in the strong coupling constant is given by 



kW) 



(3.6) 

(3.7) 
(3.8) 

The thermodynamic potential 



n f = -^{k f E f (k f )(k} - \ m }) + |m}log 



k f + E f {k } ). 



2as 

7T 



3lE f (k f )kf 



m 2 log 



k^Ejikf)^ 
' E f (k f ) 



rri4 



2fc 



3m) log 2 1 mf 



E f (k f ) 



(3.9) 



+ 6 log 



a 



Ef(kf) 



kfEf(kf)m 2 j: — m) log 



k f + E f {k f ) 



in 



f 



where 7/ = spin x color degeneracy. a s is the strong interaction coupling constant, and a is 
the renormalization scale constant which is considered as a typical chemical potential of the 
problem and chosen as ~ 300 MeV [23]. However, the EOS of hadron-quark compact stars is 
not s ensitive to the parameter a, which is numerically checked and theoretically expected from 
(|3.9p in high density region. It is suggested that the EOS of quark matter interconnected with 
hadronic matter be mainly sensitive to parameters, (B, a c ) in the current analysis. The expected 
value and restriction to (B, a c ) can be numerically extracted from properties of hadron-quark 
matter and compact stars in the sec. 6. 

The energy density and pressure are now expressed as, 



£ q = B + J2(Vf + W) 
/ 

f 



(3.10) 



where \i$ = Ef(kf) and rif are chemical potential and particle distribution for the flavor /. 
The equations of state for hadronic and quark matter, phase transition conditions and TOV 
equation [21] are employed to calculate pure-hadron, pure-quark and hadron-quark compact 
stars. 
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Masses of hadronic nuetron stars 




Central Energy Density e c (g/cm 3 ) 

Fig. la. Masses of hadronic neutron stars. Pure- 
neutron matter in linear <t-lj: M (n) max = 3.06 
Mq. Isospin asymmetric /3-equilibrium matter: 
M(n,p,e) max = 2.50 Mq. (n,p, e)-(n,p, £~, e)- 
(n,p, E~,A) matter with the coupling ratios: 
(r% N = 0.964, r% N = 1.00) and (r°_ N = 0.925, 
r%_ N = 1.00, r*_ N = 1.00) give M max = 1.98 
M Q . 
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Fig. lb. Masses of hadronic neutron stars. 
(n,p, e)-(n,p, S~, e)-(n,p, £~, A) matter with 
the coupling ratios: (r^ N = 0.677, r^ N = 2/3) 
and (r^_ N = 0.632, r«_ N = 2/3, r^_ N = 1.00) 
give M max = 1.40 Mq. 



4 The properties of hadronic and hadron-quark (H-Q) hybrid 
stars 

We have applied the conserving nonlinear a-u-p mean-field approximation with a VFC to ex- 
amine properties of hadronic neutron stars and H-Q hybrid stars. 

The properties of pure hadronic stars are produced by employing hadronic equations of 
state, phase transition conditions, self-consistent effective masses and coupling constants, and 
the TOV equation. The masses v.s. central densities of hadronic stars produced by EOSs for 
(n,p, e)-(n,p, E~, e)-(n,p, E~, A) matter are shown in Figs, la (r^ N = 1.0, r^ N = 1.0) and lb 
( r AN = 2/3, = 2/3). The solid line shows (n,p, e) matter which produces the maximum mass 
of hadronic neutron stars, M max = 2.50 M Q [18 1. LHA(n) (dot-dashed line) is the result of the 
linear o-uj mean-field approximation [13J (M max = 3.06 Mq). The arrows respectively indicate 
phase transitions from (n,p, e)-(n,p, E~, e) and (n,p, S - ,e)-(n,p, S _ ,A) matter. The EOS of 
hadronic matter with r^ N = 1.0, r£ N = 1.0 generates the maximum mass M max (n,p, A) = 
1.98 M Q (Fig. la). However, the EOS with r% N = 2/3, r% N = 2/3 as suggested by the SU(6) 
quark model, produces the small maximum mass of neutron stars, M max (n,p, S~, A) = 1.33 
Mq (Fig. lb). This is unable to support the observed masses of neutron stars. The discrepancy 
calculated by r^ N = 1.0 and r°^ N = 2/3 is clearly recognized by comparing the EOS and mass 
of neutron stars. 

All curves of £ c — M s t ar hi Fig. la and lb have simple positive slopes, dM/d£ c > 0, and one 
maximum point which is typical for lines of one critical (inflection) point, d 2 M/d£^ < 0. In 
this case, stability can be examined clearly for these simple saturating curves by the condition 
dM s tar I d£ c > [231 G3], respectively. Hence, values of positive slopes of £ c — M star are solutions 
of stable neutron stars, but the decreasing curve, dM s t ar /d£ c > 0, after the maximum point 
indicates unstable neutron stars. These characters of £ c — M star are important to examine 
stability of hadron-quark stars. 

The EOS of pure-quark matter is used to calculate masses of pure-quark stars, as shown in 
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Fig. 2. Masses of pure-quark stars with B = 80,120,150 MeV/fm 3 , (a c = 0). The stable 
pure-quark stars with M max > 1.0 are in 10 15 ~ £ c ~ 5.0 x 10 15 g/cm 3 . 
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Fig. 3. Hadron stars with VFC and pure-quark stars. The stable hadronic stars are in 5.0 x 10 
£ c ~ 10 15 g/cm 3 . 



ii 



Fig. 2. The bag constant should be B > 80 MeV/fm 3 , as normal nuclear matter becomes quark 
matter in the case of B ~ 80 MeV/fm 3 . The maximum masses of quark stars are M max = 1.59 
M Q (B = 80 MeV/fm 3 ), M max = 1.31 M Q (B = 120 MeV/fm 3 ) and M max = 1.18 M (B = 150 
MeV/fm 3 ). The £ c — M sta r curves in Fig. 2 and 3 for the current pure hadron and quark stars 
are smoothly increasing and saturating. These curves are classified as those of one critical 
(infrection) point. Therefore, stability of neutron stars in pure-hadron and pure-quark matter 
can be checked by the condition: dM stg , r /d£ c > 0, though one needs to check other conditions 
if £ c — M star curves exhibit complicated phase transitions so that critical points are larger than 

2 M- 

Stable quark stars with M s t ar ^1-0 M Q are limited in the range of central energy densities: 
10 15 < £ c < 5.0 x 10 15 g/cm 3 , as shown in bold, solid lines in Fig. 2. Stable hadronic neutron 
stars with the VFC and pure quark stars are shown in Fig. 4, to remark on the energy-density 
regions of stable stars. The results of hadronic and quark stars in Fig. 3 suggest that H-Q hybrid 
stars are possible about the central energy density, £ c ~ 10 15 g/cm 3 , as the EOS of quark matter 
could be energetically preferable at high densities, compared with that of hadrons. However, the 
overlap in central energy density between stable hadronic stars (M^ ax > 1.0) and pure-quark 
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Fig. 4. Hadron-quark stars based on isospin asymmetric (n,p, e)+VFC matter. Note that the 
quark-core of H-Q stars is stable. 

stars (Mmax ^ 1-0) is very narrow. If stable central energy densities are the same duri ng t he 
H-Q phase transition, this indicates that H-Q stars would become immediately unstable [47J. 

We have assumed that the H-Q phase transition is a first-order, and solved phase transition 
conditions for chemical potential and pressure (fin = and Ph = Pq) by employing the 
double-tangent method numerically. The TOV equation, as well as hadronic and quark EOSs 
with phase transition conditions, are applied to calculate the properties of neutron stars. First, 
we have examined H-Q hybrid stars in case of (n,p, e) + VFC and quark matter. The results are 
shown in Fig. 4. In this figure, the mass of H-Q stars is shown with a solid line and the stable 
quark core is indicated by a dotted line. The results are not sensitive to the parameter a, which 
is checked by changing a = 200 ~ 800 MeV. However, the H-Q stars are sensitive to the values of 
B and a c , B = 100 ~ 150 MeV/fm 3 , a c = 0.2 ~ 0.1 to produce the observed masses of neutron 
stars (M star ^ 1-30 M Q ). This is because the quark-EOS shifts to high energy densities if B is 
increased, meaning that the EOS becomes softer in terms of pressure, resulting in compact stars 
with a smaller mass and radius at high densities. In addition, if the QCD coupling constant a c 
is increased as a c = 0.1 — > 0.2, the masses of stable H-Q stars become small and shifted to high 
densities. Hence, appropriate values of B and a c to produce the observed data of neutron stars 
in the current EOS are correlated to each other so that if B is increased, a c should be decreased, 
such as in (B ~ 100 MeV/fm 3 , a c ~ 0.2) and (B ~ 150 MeV/fm 3 , a c ~ 0.1). The property 
of (B,a c ) in the analysis of H-Q infinite matter agrees with the results of bag-model fits to 
light-hadron spectra and renormalization group analyses in the paper by Farhi and Jaffe [49J. 

The stability of pure hadronic and quark matter is examined, respectively, by the condition 
dM star /d£ c > [23l [5l]. One should note that the stability criterion, dM star jdE c > 0, is for 
single phase compact stars, such as pure hadronic and quark stars. As H-Q stars are 2-phase 
compact stars (i.e., the quark phase for a star's core and hadron phase for a mantle), the stability 
of H-Q stars should be reconsidered. 

The H-Q stars in Fig. 4 show that the stable hadronic stars {dMn / d£ c > 0) in the central 
energy range, f c = 10 14 ~ 10 15 g/cm 3 , will undergo a phase transition (dotted horizontal line), 
reaching H-Q stars (B = 150 MeV/fm 3 ). The Fig. 4 indicates that the total mass of H-Q stars 
decreases, but a stable quark-core {dMq/d£ c > 0) develops. Therefore, this suggests that the 
H-Q star is stable, although the total mass of the star becomes smaller. Moreover, by comparing 
stable energy densities of the quark phase in Fig. 3 with those of H-Q stars, the central energy 
density of stable H-Q stars is found to be more extended for higher densities than that of single 
phase stars. This suggests that compact stars consisting of a mantle and a high density core 
are more stable than stars in a homogeneous single phase structure |49t 155] . When the QCD 
coupling constant, a c , increases, the H-Q phase transition density and quark-core will shift to 
higher densities, but the quark-core is stable and extends to higher densities. If the bag constant 
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Fig. 5a. Hadron-quark stars based on (n,p, e)- Fig. 5b. Hadron-quark stars based on (n,p, e)- 

(n,p, E~, e)-(n,p, E~, A) + VFC matter. The (n,p, £~, e)-(n,p, £~, A) + VFC matter. The 

coupling ratios: = 0.964, = 1.00) and coupling ratios: (r^ N = 0.677, r% N = 2/3) and 

(r°_ N = 0.925, r%_ N = 1.00, r^_ N = 1.00). (rg.^ = 0.632, r%_ N = 2/3, = 1.00). 

is small, such as B ~ 100 MeV/fm 3 , one will obtain a saturation curve (an inflection point for 
stability) within 10 15 ~ 10 16 g/cm 3 , as in Fig. 3. 

The H-Q stars calculated with the EOS for (n,p, e)-(n,p, S~, e)-(n,p, S~, A) + VFC to quark 
matter (B = 150 MeV/fm 3 ) are shown in Fig. 5a (r^ N = 1.0) and 5b (r^ N = 2/3), respectively. 
The maximum masses in Fig. 5a are M max = 1.61 M {a c = 0.0) and M max = 1.77 M 
(a c = 0.1); whereas, in Fig. 5b, they are M max = 1.37 M Q (a c = 0.0) and M max = 1.48 M 
(a c = 0.1). A remarkable feature of H-Q stars in Fig. 5a {r^ N = 1.0) is that they can reasonably 
explain central energy densities and maximum mass configurations [25] for stable compact stars, 
ranging from 10 14 ~ 10 16 g/cm 3 . 

In the case of the softer hadronic EOS in Fig. 5b {r l fj N = 2/3) with B = 100 ~ 150 and a c - 
correction, phase transitions occur in relatively low densities. This results in smaller maximum 
masses of H-Q stars, M max < 1.4 M . Thus, the EOS with the coupling constants (r^ N = 2/3) 
may not be appropriate to explain the observed masses of neutron stars. A physical reason from 
the hadronic sector is clear in terms of the EOS; however, it is very interesting to investigate 
how hadronic and quark models reconcile the problems pointed out in the paper. 

If the QCD coupling constant is set large (a c > 0.1) in the soft EOS of Fig. 5b {r l f {N = 2/3), 
generation of the quark-phase will move up to a very high density, separating the energy densities 
of hadron and quark phases completely. In this case, no stable hadron-phase can exist or coexist 
even in a crust of the surface of quark stars, resulting in literally pure-quark stars. Although 
more precise and detailed analyses are needed for many-body interactions of hadronic and quark 
matter, the coupling constants for hadrons (r^ N = 2/3 or r^ N = 1/3), as suggested by the SU(6) 
quark model, would not be appropriate to examine the properties of nuclear and hyperonic 
matter, or the maximum masses of high density compact stars. 

5 Concluding remarks 

The quantum hadrodynamics (QHD) and mean-field approximations are interesting for de- 
scribing hadronic many-body systems at large distances, but they must ultimately break down 
at short distances where QCD is valid; hence, it is stimulating to investigate how properties of 
hadron and quark dynamics will interconnect and influence each other dynamically [13\ I14j . We 
have applied the conserving nonlinear a-uo-p mean-field approximation to study self-consistent, 
density-dependent interactions of hadrons and examined hadronic dynamics by employing cou- 



12 



pling constants required from the QCD degrees of freedom. Discrepancies between hadronic 
and QCD predictions are shown in terms of effective masses, effective coupling constants, in- 
compressibility and symmetry energy, saturation properties of hyperon binding energies and 
maximum masses of high density compact stars. One of our purposes is to compare predictions 
indicated by effective hadronic and quark models so as to clarify interrelations and distinctions 
between them. 

As saturation properties of symmetric nuclear matter are self-consistently related to those 
of hyperons, reproducing saturation properties of hyperons should be one of fundamental prob- 
lems of nuclear many-body approximations pQ — [7]. Although more quantitative and model- 
independent analyses are needed, it is concluded that hyperon coupling ratios required by SU(6) 
quark model for the vector coupling constants [37\ I38j are not appropriate in order to generate 
saturation properties of hyperons. The coupling ratios of hyperons, r°^ N ~ 1.0, are appropri- 
ate to explain Fermi-liquid properties, saturation of hyperons [181 I19j and maximum masses of 
neutron stars. Discrepancies of hyperon coupling ratios required by hadronic and SU(6) quark 
models should indicate that both effective approaches be improved with each other as consistent 
theories for hadronic physics. 

Heavy-ion collision experiments as well as neutron stars are useful to examine applications 
of quark models and constraints for hadronic calculations |56} 157] . since conditions of hadron- 
quark phase transition depend on both equations of state for hadronic and quark matter. The 
appropriate values of (£>, a c ) can be independently suggested from properties of infinite matter, 
which will qualitatively support properties discussed in the paper |49j . The stable pure-hadronic 
stars and pure-quark stars based on MIT-bag model exist in different energy densities respec- 
tively, but if a first-order phase transition is assumed, the energy density of stable hadron-quark 
stars expands from 10 14 to 10 16 g/cm 3 . In the range of energy density, pure-hadron stars are less 

than M(n,p, e) max < 2.5, but H-Q stars are Mmax > 1-0 M & . It reasonably explains expected 
masses and energy density relations, M — £ c |25j . Therefore, existing high density compact stars 
are more likely to be regarded as hadron-quark stars rather than pure-hadron and pure-quark 
stars. These results should be further examined quantitatively from empirical data and effective 
theories of hadrons. 

It is noteworthy that density-dependent many-body effects simulated by nonlinear a-u-p 
interactions are more important than those of VFC in the nonlinear a-u-p mean-field approxi- 
mation. The contributions of VFC should be investigated further in more complicated nonlinear 
a-oj-p HF, BHF approximations in order to study vacuum-fluctuation and density-dependent 
corrections. The role of chiral symmetry in hadronic models should be examined by extending 
the current nonlinear a-uo-p mean- field approximation by including it- meson, which would reveal 
physical meaning and significance of chirality for hadronic and quark models. 
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